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Abstract 

Employing standard results from spectral geometry, we provide strong evidence that in the clas- 
sical limit the ground state of three-dimensional causal dynamical triangulations is de Sitter 
spacetime. This result is obtained by measuring the expectation value of the spectral dimension 
on the ensemble of geometries defined by these models, and comparing its large scale behaviour 
to that of a sphere (Euclidean de Sitter). From the same measurement we are also able to con- 
firm the phenomenon of dynamical dimensional reduction observed in this and other approaches 
to quantum gravity - the first time this has been done for three-dimensional causal dynamical 
triangulations. In this case, the value for the short-scale limit of the spectral dimension that we 
find is approximately 2. We comment on the relevance of these results for the comparison to 
asymptotic safety and Hofava-Lifshitz gravity, among other approaches to quantum gravity. 



1 Introduction 

The path integral approach to the non-perturbative quantisation of gravity has made unprece- 
dented progress in recent years. While a Euclidean path integral approach was suggested and 
studied many years ago [T], most of the non-perturbative results have remained at a formal level, 
because of the ill-defined nature of the functional integral. Discretisation offers a way to better 
define such a path integral, as a continuum limit of a "lattice" theory, in analogy to standard lat- 
tice quantum field theory. The causal dynamical triangulation (CDT) approach to discretisation 
[2j [3] has proven to be very useful in this regard. In particular, it has allowed the calculation of 
expectation values of physically relevant observables as continuum limits of discrete observables. 
These results provide some evidence that the CDT theory is approximating general relativity at 
large scales, the first step to confirming that a theory of quantum gravity has been discovered. 

In particular, several results are consistent with the conjecture that the CDT theory is able to 
produce average geometrical configurations that approximate to the de Sitter solution of general 
relativity at large scales. This is consistent with the claim that CDT theory has a semi-classical 
regime with the expected behaviour. Not long ago, results on four-dimensional CDT [5] showed 
that the effective large-scale dimension of the universe was approximately four, a non-trivial result 
in view of the possibility of pathological geometries dominating the path integral (as in the previ- 
ous Euclidean DT approach [5]). Measuring the expectation value of the volume of spatial slices 



as a function of cosmological time shows a good match to de Sitter geometry [6]. Fluctuations 
about de Sitter space in the CDT model can also be matched by a reduced "minisuperspace" 
quantum cosmology model derived from general relativity. Such results give the first evidence 
that, not only is the theory producing a solution to general relativity, but also the dynamics is 
consistent with what we know of quantised general relativity. 

Still, it is difficult to measure a large and general class of geometrical observables with the 
computational methods that are being applied. The results confirming de Sitter space are so far 
limited to information on the size of spatial slices, and dimension. The difficulty is in finding 
genuine geometric observables (in keeping with the symmetries of general relativity) which are 
relevant at large scales (in keeping with the coarse-graining that we understand to be necessary to 
obtain physical results in quantum field theories). One approach to this problem was the proposal 
in [7] of a coarse- graining scheme for CDT, but the application of the method to simulations has 
not yet been achieved. It would be desirable to identify a set of observables that could be 
more easily measured in Monte Carlo simulations of the CDT model, and rich enough to fully 
characterise the geometry. This would be necessary in order to have a complete test that the 
model can produce an approximation to de Sitter at large scales (and, very similarly, that the 
geometry in the classical limit Lpi anc k — > is exactly de Sitter). 

One of the dimension estimators, the "spectral dimension", gives a hint of how this might be 
done. In [S] it is explained how a scaling effective dimension estimator can be derived from the 
spectrum of the Laplace-Beltrami operator on a geometry. The spectral dimension is a function 
of a "diffusion time" parameter corresponding to the relevant scale, and the function is in 1- 
to-1 correspondence with the spectrum of the Laplacian. Its values can be found by examining 
diffusion processes on the geometry, which can easily be discretised on a CDT. Interestingly, 
the results that can be gained in this way go beyond the dimension; some geometries (such as 
the Euclideanised version of de Sitter, the sphere) are entirely characterised by this spectral 
information. Thus, by a closer examination of this information, readily available from Monte 
Carlo simulations, it may be possible to find strong evidence for the emergence of de Sitter space. 
Indeed, this could be as close to a necessary and sufficient test as it is possible to get, within the 
practical limitations of computer simulations. 

As well as the large-scale considerations, it is interesting to consider results from spectral 
geometry at short scales. In 4D, current results [8] suggest a kind of "dynamical dimensional 
reduction", in which the spectral dimension approaches 2 at short scales. This gives an example 
of the highly non-trivial nature of the UV regime of quantum gravity, gives some hints of how 
gravity might be a consistent theory in the UV limit, and also provides an observable to compare 
to results from other approaches to quantum gravity. 

These ideas are explained in more detail below. Applying them to the 3D CDT model, we 
find that large scale spectral dimension function derived from simulations matches the expected 
classical function very well at large scales. We also find that the discrepancy is reduced as we 
approach a continuum (classical) limit. As well as this, a new dimension estimator, based on the 
scaling of the spectral dimension function with the number of simplices N in the CDT, gives a 
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result consistent with three, with good accuracy. 

These ideas are applied to the 3D model of CDTs here, for two reasons. The first is that the 
simpler model allows simulations of larger linear size than the 4D model, and provides a testing 
ground for the ideas, which seem to be equally applicable to the 4D case. Secondly, measurements 
of the spectral dimension have not previously been made in the 3D case. This is of interest in 
itself, as measurements of the small scale spectral dimension for models of different dimensionality 
may allow a better comparison between the CDT models and various other candidate theories 
for quantum gravity, some of which have predictions for the small scale spectral dimension. Here, 
we find the small scale dimension of 3D CDTs to be approximately two, as for the 4D model. 

After briefly reviewing the essential concepts of spectral geometry in Sec. [21 and the basics of 
CDT in Sec. O we describe in Sec. H]how the two are combined. Finally in Sec. Owe present our 
results and discuss their consequences. After the Conclusions, two appendices are added with 
some supporting material. Our results lead to a strong confirmation of the de Sitter character of 
the large scale geometry and some interesting comparisons of quantum gravity theories in three 
dimensions. 



2 The heat trace and geometry 

First, we review some relevant results from spectral geometry. Consider a d-dimensional closed 
Riemannian manifold M with a smooth, fixed metric g^ u . The heat kernel K(£q, £, a) is a function 
on M x M x M_|_, which solves the heat equation, 

j^M^o, a) + A g K g (to, £,a) = Q, (2.1) 

with initial condition 

if g (^,a = 0) = ^ ( lzff . (2.2) 

Here, a is the "diffusion time" (not to be confused with any physical time when considering 
the manifold as a spacetime), and A g = — p^V^Vj, is the scalar Laplacian on M (acting on £), 
where is the covariant derivative compatible with gw. In mathematical terms, the heat kernel 
K g (^Q,^, a) is the Green function of the heat equation, while its physical interpretation is as the 
probability density of diffusion from £o to £ in diffusion time a. The solution of (|2.ip is formally 
given by K = (^\e~ aA \£,o) , or in terms of eigenvalues Xj and eigenfunctions </>j(£) of A g 

K g fa,M = xy^-ceweo) > ( 2 - 3 ) 

3 

where it has to be understood that in the case that the spectrum is continuous the sum would 
be replaced by an integral. 

From this, we define a function that is a generally covariant observable, the heat trace, 

W = ^[ ^^MY\K g ^^,a), (2.4) 
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where V g = J M d£ y/\g(£)\ is the volume of the manifold, a normalization factor that is introduced 
for convenience. The heat trace is the probability density for returning to the starting point of 
diffusion at time a, integrated over M. It is determined by (and determines) the spectrum via 

w4^ e " v ' (2 - 5) 

9 3 

where the eigenvalues of A g are repeated according to their degeneracy. From this we see that 
the value of the heat trace at diffusion time a is not significantly affected by eigenvalues much 
larger than 1/a. This is the relationship between diffusion time and the scale being "probed". 
The heat trace is also related to invariants of the geometry through the well-known heat-trace 
expansion [9J, 

1 oo 

where a n is a series of curvature invariants, the first three of which are 

ao= [ dtVWfi, Oi = i / d£ V\g~W\R(0 (2-7) 

JM O JM 

« 2 = I d ZVW)\{$R 2 ~ 2R^R^ U + 2R fMUpT R^} , (2.8) 
ooll J M 

where R, R^v, and R^ U pr are the scalar curvature, the Ricci and the Riemann tensors respectively. 
Thus, for these smooth geometries, the heat trace determines dimension, volume, average scalar 
curvature and other curvature invariants. In fact, in the case of the 3D and 4D spheres, the 
heat trace completely determines the geometry: it is known that if M is a closed, connected 
Riemannian manifold of dimension 2 < n < 6 with smooth metric g^vi and if (M,g) has the same 
spectrum as the n-sphere S n with the standard metric, then (M,g) is in fact isometric to S n [TP] . 
This is not true for all geometries, however. If two geometries have the same heat trace function, 
they are called isospectral. 
The function 

can be used to find the dimension of the manifold. Considering for example a flat space, the heat 
trace reduces to 

P g (a) = (47Rx)- rf / 2 , (2.10) 
and so we see that d s (a) = d in this case. In the general case, 

d s (o) = d-2 ^Z inana n , (2.H) 

and we see that the spectral dimension reduces to the topological dimension only in the limit 
a — > 0, with a slope determined by the total curvature. At very large diffusion time, one can 
see from (|2.5|) that the spectral dimension has an exponential falloff determined by the lowest 
eigenvalues. 
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Up to here we have recalled well-known facts about spectral properties of classical manifolds. 
In the context of quantum gravity we expect that at very short scales spacetime will not look 
classical. It will instead be replaced by some structure which, although unknown, has for some 
time generally been known as "spacetime foam". In order to make some sense of this general 
expectation, one has to study in detail some specific model of quantum gravity and try to find at 
least an effective description of what happens to spacetime when we probe it at shorter and shorter 
scales. In this vein, we need to find geometrical "probes" from which to obtain some indications 
of what is going on at the Planck scale. One possibility is to take, as such a probe, the heat 
kernel trace (adapted to the quantum gravity model in consideration) and study its deviations 
from classicality. This way we can gain some interesting insight on the short scale effects and 
also test whether the model has good classical properties at large scales. In particular we can 
use (|2.9p to define an effective notion of dimension for the quantum geometry. In this case, we 
expect that the spectral dimension function of a good quantum gravity theory will approximate 
well to that of a classical geometry at large diffusion times, while at small diffusion times there 
will be significant deviations. Below, we match the average spectral dimension function we derive 
from CDT simulations to that of a candidate classical geometry at large diffusion times, and then 
examine the small diffusion time behaviour to find properties of the spacetime foam. 

Some of these ideas have been applied to d = 4 CDT models by Ambj0rn et al. in [51 [3], using 
methods that we compare to our own in section 15.31 The results they found showed for the first 
time a scale-dependent effective dimension which agrees with the classical (topological) one at 
large scales but reduces to 2 at short scales. This dimensional reduction is suggestive of a picture 
of quantum gravity in which the theory self-regularizes its behaviour, making itself safe from the 
UV-catastrophe which is usually associated with quantum gravity. Indeed d = 2 is the critical 
dimension of gravity, and the theory presents no problems in that case. Such a picture seems 
compatible with that of the "asymptotic safety scenario" [12"| [13] , for which the existence of 
a non-trivial fixed-point of the renormalization group equations implies that the theory near the 
fixed-point behaves like a two-dimensional theory (as deduced by the anomalous dimension). The 
spectral dimension has also been derived in the asymptotic safety scenario |14j and the results 
agree with CDT for the d = 4 case. 

Finally the spectral dimension has been studied also in the context of spaces with quantum 
group symmetry [T5] , in loop quantum gravity and spin foams W7\ [TB] , in Hof ava-Lifschitz 
gravity |19j . and in the strong-coupling limit of the Wheeler-DeWitt equation [20 1 . We will 
comment more on the relation among these results, and their relation to our new findings, in 
section 15.41 

3 A brief survey of CDT 

The CDT models are a concrete proposal to define a path integral for gravity. A motivation 
for seeking such a definition is of course the hope that gravity might make sense if defined non- 
perturbatively. The traditional way to pursue a non-perturbative evaluation of a path integral in 
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quantum field theory is to replace the continuum spacetime with a fixed lattice, a procedure which 
allows one to do actual computations (typically numerical ones, via simulations) with an object, 
the path integral, which otherwise has no meaning. Then one tries to recover the continuum 
theory, heavily borrowing concepts and procedures from the theory of critical phenomena. 

The CDT approach exactly parallels the usual lattice field theory with one fundamental 
difference: the fixed lattice is replaced by an ensemble of random triangulations. This is required 
by the fact that gravity is a theory of dynamical geometry, with no background spacetime fixed 
a priori. 

More concretely, one defines an ensemble of "triangulations" to work with, a triangulation 
being defined by a simplicial manifold, i.e. a collection of d-dimensional flat simplices (the gen- 
eralization of triangles and tetrahedra) glued along their (d — l)-dimensional faces and such that 
the neighbourhood of any vertex is homeomorphic to a <i-dimensional ball. A dynamical tri- 
angulation is one in which all the simplices are taken to be equilateral. In the simulations we 
usually work with dynamical triangulations having a fixed number of d-simplices N, which we 
will denote T/v- The ensemble of such triangulations {Tjv} is obtained by gluing the iV simplices 
in all possible ways allowed by the simplicial manifold condition. Furthermore, to avoid the sick 
behaviour that was found in the old models of dynamical triangulations, CDT models have one 
further restriction on the ensemble: only triangulations with a global time foliation, with respect 
to which no topology change occurs, are allowed. For more details on the geometrical meaning 
of this restriction and on its implementation see |21j . 

Once the ensemble is specified one can construct the partition function (Euclidean version of 
the path integral) as 



where S(Tn) is the bare action, and C(T/v) is the order of the automorphism group of Tjv, a 
symmetry factor naturally appearing when summing over unlabeled triangulations. Since we 
wish to recover general relativity in the classical limit, it is customary to use as a bare action the 
Einstein-Hilbert action adapted to a simplicial manifold, which is known as the Regge action. On 
a dynamical triangulation, the Regge action takes the very simple form 



where Kd and Kd-2 are two coupling constants depending on the cosmological and Newton's con- 
stant appearing in the Regge action, and Nd—2 is the number of (d — 2)-dimensional subsimplices 
(also called bones or hinges). 

In principle one could use a different action, with more parameters, but at this stage this would 
only complicate the analysis of the results, and in a minimalist attitude such a generalization of 
the CDT models is usually postponed till the moment (if ever) at which the model itself will ask 
for such an extension. For example in 3 + 1 dimensions a new parameter has been introduced in 
the action, without which no physically interesting region would exist in the phase diagram [3]. 
Furthermore, we need to remember that, as a consequence of topological relations, only d/2 (for 




(3.1) 



N T N 



S{T N ) = K d N - Kd- 2 Nd-2 



(3.2) 
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d even) or (d + l)/2 (for d odd) among the values {No, N\, ...Nd-i, N} are independent. Hence, 
for d = 3 and d = 4 only two of such variables are independent, and as a consequence, if we 
want to keep the action linear in Nj, we only have two coupling constants. The counting changes 
if some anisotropy is introduced in the model, by assuming that the ratio a = If/lg between 
the lengths of time-like and space-like edges^ is different from one. In such case, one finds that 
additional variables are needed in order to keep track of the orientation of the subsimplices, and 
new topological relations are found too. The counting for the anisotropic models was carried out 
in [21j, and one has that for d = 4 there are 10 variables and 7 constraints, leaving 3 independent 
variables, a fact that was used in [3] to introduce the new parameter. In d = 3 the situation is 
instead unchanged with respect to the isotropic there are 5 constraints for 7 variables, 

and hence again only 2 independent variables. For this reason it does not make sense to introduce 
in 3D the analogue of the new parameter used in 4D. 

In this paper we concentrate on the case d = 3, for which very few analytical results are 
known [22\ [23] because of the difficulty in solving statistical models in dimensions higher than 
two. Hence we will resort to the method of Monte Carlo simulations. In the simulations we will 
use the topological constraints to trade the variable N\ for Nq, which is easier to keep track of, 
and replace (|3.2|) for d = 3 by 

S(T N ) = k 3 N - k N . (3.3) 

Furthermore, as we mentioned, in the computer simulations we work at fixed volume, and hence 
we replace (|3.ip by 

where we have made use of the simple form of the action (|3.3j) . Note that the partition function Z 
is the discrete Laplace transform of Zjy with respect to N. The expectation value of an observable 
A is calculated as 

T N 

which is related to the expectation value as a function of K3 via 

(A) = ^e'^ N Z N (A) N . (3.6) 

N 

Note that all the quantities appearing in (|3. 2|) are dimensionless. Dimensions can be reintro- 
duced in terms of the edge length a of the simplices, the equivalent of the standard lattice unit, 
which is necessary when talking about the continuum limit. An essential part of the continuum 
limit procedure is to have A" — > 00 and a — > in such a way that the physical volume 

V ~ a d N (3.7) 

We are using here a Lorentzian language even though our signature is Euclidean; this is possible as in CDT 
it is always clear which edges in the triangulations are to be thought as time-like edges, and a Wick rotation is 
possible, at least before the continuum limit is taken. 
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remains finite. This implies that when we want to give dimension to a dimensionless quantity by 
multiplying it by a n , in practice we multiply it by N~ n / d . In the continuum limit we expect that 
large-scale observables will become independent of the cutoff a, hence we expect to see finite size 
scaling when working with simulations at sufficiently large N, i.e. we expect that an observable 
/(<7i, ■■■,g m ; N), depending on a set of m couplings or variables {gi}, will satisfy 

f(g l N-li,...,g m N-^;N) = f(g 1 N'-ir,...,g m N , -^;N') (3.8) 

The natural expectation would be that (|3.T|) and (|3,8p hold with {n^} given by the expected 
length-dimension of the {gi}, but this is not guaranteed a priori and it is instead used as a check 
of the good classical properties of the model, as we do in the following. 

4 Application of spectral geometry to CDTs 

The expectation value of the spectral dimension function for CDTs was found from Monte Carlo 
simulations. Finite size scaling analysis is applied to the results, as in |24j . Some previously 
existing code for the Monte Carlo simulations (used in [24J) was adapted for this purpose. These 
simulations fixed the spacetime topology to S 2 x S 1 , i.e. spherical spatial sections and cyclical 
time. The cyclical time is merely a convenience; it has been seen that the simulations produce 
an extended "universe" part consistent with a 3-sphere, and a "stem" part in which spatial slices 
remain near the minimum set by the lattice scale. The "stem" part is thought to be a discreteness 
artifact caused by the topology used in the simulation. Values of N up to a maximum of 200k 
(meaning 2 x 10 5 ) were studied, although some errors for the larger values of N are greater since 
less configurations could be generated to be averaged over, within practical time constraints. 

All simulations were carried out with coupling constant kq = 5, in the "extended phase" of the 
CDT phase diagram, where previous evidence points to the emergence of well-behaved geometry. 
The total number of time-steps was set to T = 96. We approach the continuum limit with ko 
fixed, N — > oo and keeping the dimensionful volume V fixed. In this case the Planck length is 
fixed in terms of lattice units, and as a result the dimensionful Planck length must go to in this 
limit, and so this represents a classical limit. In order to keep the Planck length finite a different 
limit would have to be taken with kq a function of N [2]. 

The heat trace of a CDT configuration can be found by appropriately discretising the diffusion 
process. This process uses a probability distribution defined over the set of simplices of the CDT. 
The analog of the continuum diffusion equation (|2.1|) is 

Kt{x,x q ,(t + 1) = (1-x)Kt{x ) xq ) (t) + - ^2 K T (x',x ,a) (4.1) 

x'eG(x) 

where x and x' labels of simplices, G{x) is the set of simplices glued to x (of which there are 
4 in the 3D case), and x plays the role of a diffusion constant, which can be between 1 and 0. 
The distribution is initially peaked at Kt(x,xq,0) = 5 x ,x - Similarly to the continuum case, for 
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a particular triangulation T the average probability of return is 



Pt{t) 



1 



(4.2) 



N(T) 



and the analog of the spectral dimension of eqn. (|2,9p . d s (cr,T), is defined by replacing the dif- 
ferential with a finite difference in the obvious way. The Monte Carlo simulations estimate the 
ensemble average of this quantity, 



The time normalisation makes no difference to our results at this stage. Similarly the value of 
X is set purely from practical considerations, to minimise discreteness artifacts while allowing 
diffusion to proceed fast enough to probe large scales within the possible number of diffusion 
steps. As explained in [3], at very early times the value of Pt has an oscillatory behaviour due 
to discreteness effects. This is somewhat mitigated by reducing the value of x> which was chosen 
to be 0.8 for the simulations. It is also possible to introduce a different rate of diffusion through 
spacelike and timelike triangles on which the simplices are glued, which amounts to a scaling of 
timelike and spacelike distances, but this was not done here for reasons explained below. 

In principle D s would be calculated from (|4.3p by carrying out the diffusion process for all 
tetrahedra in each Monte Carlo configuration. However, a good estimate of this result can be 
found by randomly sampling the tetrahedra x at which to calculate the probability of return 
Kt{x, x, a). One tetrahedron was used as starting point for each configuration, selected uniformly 
at random from all tetrahedra not in the stem, and this was found to give reasonably small errors. 

Before the process was simulated on the CDT, the CDT structure was converted from the 
data structure optimised for Monte Carlo moves to one optimised for the diffusion, resulting in a 
considerable speed-up. This allowed the diffusion process to be followed for of the order of 20 000 
steps for all configurations, in the order of hours of computer time, more than enough for our 
purposes. The diffusion process was carried out with the stem excluded. The stem here is defined 
as including all simplices in "small" slices (< 20 simplices), and slices in ranges that lie between 
two small slices, but do not include the maximal slice (which would indicate that this range is the 
universe part). The value of N used in scaling calculations should, strictly speaking, be altered 
to take this into account, but this was not found to alter results in most cases q 

One of our aims here is to compare these results from CDT simulations to the de Sitter 
geometry that we expect, from previous studies [23], to find at large scales. The first step, then, 
is to determine the spectral dimension function for de Sitter, or the Euclideanised version, the 
sphere. However, in practice, there is a complication here: the lattice structure of CDTs breaks 
the symmetry between space and time. Therefore we should not expect the ratio between the 
lattice lengths in the space and time direction that is used to define the bare action (set to 1 in 

2 The exception is the value of the scaling dimension given in equation (|5.2p . There, the value given was calculated 
using the total number of simplices less the average number of simplices in the stem, which was found to alter the 
result slightly. Elsewhere the values of N used are the total numbers of simplices. 



D s (a,N) = {d s (a)) N = — 




d s (a,T N ) . 



(4.3) 
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this case) to be the correct one to use when comparing to a continuum geometry like de Sitter. 
We are at liberty to perform some global rescaling of the time t = ai, where % is the discrete time 
step, as was done in [2] for the case of time-dependence of the volume of the spatial slices in 4D. In 
that case, the time t was taken as a coordinate time so that the scaling to the "true" cosmological 
proper time could be expressed as r = yfgut. This becomes important when considering how to 
compare to de Sitter. 

In the case that the observable under investigation is the spectral dimension, for which we 
don't know the explicit dependence on the proper time, there are two ways to proceed. The 
first, roughly, is to scale the CDT before comparing to the sphere. This means carrying out the 
diffusion process given above with different diffusion rates for spacelike and timelike faces of the 
simplices. The hope would be to extract the spectral dimension plot for a sphere by appropriately 
setting the ratio of these diffusion constants. However, there is a limit on how fast diffusion in 
the time direction can be made using this technique. The reason is that no simplicial path across 
more than two time slices in a CDT traverses only spacelike triangles, as each tetrahedron has 
at most only one spacelike face. Because of this, no matter how fast diffusion across spacelike 
triangles is, diffusion across many time slices will not become arbitrarily fast. It was found in 
simulations that the scaling made possible in this way was not sufficient to reach the sphere, and 
so only results with equal diffusion rates on timelike and spacelike edges are given belov\l§. 

The second way to proceed is to scale the sphere we are comparing to, rather than the CDT 
itself. That is, we carry out the diffusion process using equal diffusion constants on spacelike and 
timelike triangles, but scale the cosmological proper time in the sphere that we wish to compare 
to before computing its spectral dimension function. This means that we must compare to the 
spectral dimension function for the following metric, which we call the "stretched sphere": 

ds 2 s z = r 2 {s 2 d^ 2 + sin 2 ijj (d6 2 + sin 2 9 dcp 2 )) , (4.4) 

where ip,6 S [0,7r] and 4> £ [0, 27r], and s is the deformation parameter (for s = 1 we have the 
standard metric on a sphere of radius r). The spectrum for this geometry is derived in App. lAl 
and the spectral dimension can be computed from it, summing numerically the series (|A.14p . 

Our second goal will be to investigate the short-scale behaviour of the spectral dimension 
function, in order to confirm in 3D the phenomenon of dynamical dimensional reduction observed 
in 4D [HIE]. Going beyond the qualitative result, it is interesting to obtain a quantitative estimate 
of the asymptotic value (for a — > 0) of the spectral dimension in order to compare it with results 
from other approaches to quantum gravity. 

3 However, it is interesting to note that this would not be a problem if the diffusion was carried out on the vertices 
of the CDT lattice itself rather than on the tetrahedra; in that case there are paths traversing many time-slices 
consisting only of timelike edges. 
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5 Analysis of the results 



5.1 Dimension from Scaling 

Before comparing to the sphere, the first task is to verify that the expectation value of the 
spectral dimension function is scaling with N as might be expected for an extended 3D geometry. 
Examining the scaling of the diffusion time of the spectral dimension functions with N gives a way 
to estimate the large-scale dimension, similarly to the scaling dimension derived from volume- 
volume correlations discussed in [23]. Here, on large scales (i.e. large diffusion times) we expect 
the functions to scale according to 

D s (aN- 2 / d , N) = D s (aN'- 2 / d , N') , (5.1) 

because the diffusion time should scale with the square of the linear scale, from dimensional 
considerations (see eqn, (|2.ip ). We do not expect such simple scaling at small diffusion time, as 
this is the regime in which quantum effects become important, and these are probably governed 
by an additional parameter, most likely the Planck length. Since we keep kq fixed we know that 
the Planck length will remain fixed in units of the cut-off, and so we cannot apply the finite-size 
scaling reasoning of the end of Sec. [3] at these short scales. 

Figure [T] shows that a good match between data sets with different value of N is possible with 
d = 3, and illustrates that the simulations are showing the expected convergence to a continuum 
limit at large diffusion times. The best estimate of d was found by comparing the ./V = 50k and 
N = 100k data sets using best overlap methods^- The range of times to be compared was fixed 
for the N = 100k data to 2 < eriV -2 / 3 < 7.4 to exclude short scale quantum effects. The data 
was scaled appropriately for various values of d and compared using linear interpolation. The 
best overlap was found at 

d = 2.99 ±0.12, (5.2) 

where the errors are the Monte Carlo random errors (finite difference errors and so on being 
negligible compared to these). This method gives an estimate of the dimension that improves on 
the accuracy of previous methods, and gives yet more evidence for the 3 dimensional nature of 
the geometry at large scales. In view of this result, we define a = criV~ 2 / 3 as an appropriately 
scaled diffusion time, so that we can compare simulations at all values of N to a stretched sphere 
of fixed volume. 

5.2 Comparison of CDT results to the stretched sphere 

Our main goal here is to determine whether the spectral dimension function for a stretched sphere 
can match that of eqn. f|4.3j) at large diffusion times. Figure [2] shows that a scaled sphere with 
r = 1.20, s = 1.96 matches the N = 200k Monte Carlo data for a > 1.5. The good quality of 
the 2-parameter fit is evident from the plot. Note that the quality of the fit is not generic in any 

4 Comparisons to N = 200k data were also consistent with d = 3 but the errors were larger due to larger Monte 
Carlo errors. 
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a/if 13 

Figure 1: The Spectral dimension data against scaled diffusion time a — aN~ 2 / 3 . Data for N =50k, 100k 
and 200k is shown in blue, green and red respectively. Error bars are suppressed on most data points for 
clarity (this is done throughout the paper). The curves should converge to a limiting curve as N — > oo. 
At these values of N it there is good evidence of convergence to a continuum limit: the curves for N = 
100k and N = 200k agree within error for a > 0.5, as do the N — 50k and N = 100k curves. 

sense: a simple exponential gives a bad fit over this range, and altering the value of s by more 
than 0.01 also makes a noticeable difference to the quality of the fit. Evidence was also found 
that the spectral dimension function for CDTs converges to that for the scaled sphere as N — > oo, 
as shown in figure 

This remarkable fit not only shows consistency with the emergence of a 3-dimensional extended 
geometry on large scales, but also for the consistency between all the spectral properties of the 
scaled sphere and those of the results of the CDT simulations at large scales. In other words, 
these results provide strong evidence that the average CDT geometry for the largest simulations 
is isospectral to some geometry that approximates to the stretched sphere at large scales. It also 
indicates that, in the classical limit, the geometry is isospectral to the scaled sphere. We interpret 
this as the strongest evidence so far that the 3D CDT simulations are approximating 3D de Sitter 
space at large scales. 

It is interesting to note the form of the quantum correction to the classical geometry, as given 
in figure [2j At very low times the spectral dimension for the CDT data is lower than the classical 
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Figure 2: The Spectral dimension data from CDT simulations with N — 200k against scaled diffusion time 
a = aN~ 2 / 3 , plotted in black, fitted to the spectral dimension plot for the scaled sphere with r = 1.20, 
s = 1.96, superimposed in green. The two-parameter fit agrees with the data for a > 1.5. at very low 
times the "quantum correction" to the posited classical geometry is negative, but it becomes positive in 
an intermediate range. 

value, as found in the 4D case (this is investigated more fully in the next section). However, there 
is an intermediate regime in which the quantum correction is positive, and persists for a longer 
diffusion time than might be expected from looking at the behaviour at very low diffusion times. 

5.3 Results at shorter scales 

These results have implications for previous methods used to measure the large scale dimension 
of CDTs. In [H], in the 4D case, Ambj0rn, Jurkiewicz and Loll (AJL hereafter) take a short range 
of early diffusion times, and fit the CDT spectral dimension data to the function 

D s (a) = a--^—. (5.3) 
a + c 

This fit provides estimates of the dimension as a — > and as a —> oo. In particular, the estimation 
of the large scale dimension rests on two assumptions. The first is that the finite size effects have 
a negligible effect on the fit over the range fitted. The second is that the fit is good into the range 
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1/ N 



Figure 3: The overlap between the CDT data for D s (a,N) and that for the scaled sphere with r = 1.20, 
a = 1.96, as a function of 1/N, with errors. Data for N — 70k , 100k, 140k, and 200k is used here. The 
difference between the two functions is measured as the integrated absolute difference between them in 
the range < a < 7.4. As 1/N — > the difference comes closer to 0. The results are consistent with the 
convergence of the CDT data with the scaled sphere spectral dimension function as N — > oo. (At smaller- 
values of N than 70k, there is competition between the positive and negative quantum corrections that 
obscures the convergence.) 

at which the quantum corrections to the spectral dimension become negligible. In other words, 
the classical spectral dimension function being limited to as a —> oo is approximated as a constant 
function, and it is assumed that the quantum corrections affect only the non-constant term. As 
N is increased, and the range of the fit is kept as a constant range in a, the (negative) finite 
size errors become smaller. However the errors from the quantum deviations from the classical 
geometry will not go to zero unless these quantum deviations are well- fitted by b/(o + c). 

The results of the previous section show this latter assumption to be only approximately 
correct, in the 3D case at least. To illustrate this further, we repeat an analysis similar to the 
AJL analysis. At N = 70k, we find that in the 3D case the data is similarly fitted with 3 
parameters, with the large scale limit being a constant. In this case, the following fit is better 
than the rational fit given above: 



The fit is taken between a = 20 and a = 300 to avoid both short scale discreteness effects and 
finite size effects, as in [8J. As can be seen from Fig. U the fit is good (although it matches 
less well at large and small times), and the estimates derived (with errors calculated combining 
Monte Carlo random errors with fitting errors and finite difference errors) are the following: 



D s (a) = a + be 



—cx 



(5.4) 




(5.5) 
(5.6) 



D s (0, N = 70k) = 2.12 ± 0.04 . 
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• CDT data N=7Q 000 
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Figure 4: The Spectral dimension data with number of simplices N = 70k, at small diffusion times, is 
plotted with error bars in black. The exponential fit is superimposed in a lighter colour. 

The large scale dimension is consistent with 3 here, even with the assumptions given above, 
and one expects the finite size error to be small and negative. This is similar to the situation 
found by AJL in the 4D case, where the large scale dimension is consistent with 4. However, as 
one goes to yet higher N, the technique begins to slightly overestimate the dimension in the 3D 
case. 

At N = 200k we have a similarly good fit with the same function (see Fig. [5]). The estimates 
in this case are 

D s (oo,N = 200k) = 3.05 ±0.04 ; (5.7) 
D s (0,N = 200k) = 2.04 ±0.10 . (5.8) 

Here, the value 3 for the large scale dimension begins to looks doubtful. This is explained 
by looking at the spectral dimension plot at larger scales, as in figure [5J There we see that the 
quantum correction to the spectral dimension function for the conjectured classical geometry 
actually becomes positive, and stays significant for some time after the maximum value of D s is 
reached. But the above fits are taken over a short range before the maximum is reached. This 
means that the assumption that the fit will have a large scale limit free of quantum effects is 
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• CDT data N-200 000 



50 100 150 200 250 300 



Figure 5: The Spectral dimension data with number of simplices N — 200fc, at small diffusion times, is 
plotted with error bars in black. The exponential fit is superimposed in a lighter colour. 

not exact, giving a positive bias in D s (oo). However, the correction to the estimates given above 
because of this effect cannot be larger than the maximum quantum correction to the classical 
spectral dimension, which is 0.28 . Thus, we see that the extra error that we should include is not 
so large that it totally destroys usefulness of the technique; even with the approximate assumption 
this method does show that the large-scale dimension converging to a value consistent with 3, with 
approximately 10% error. The dimension estimations of the previous section also give accurate 
new confirmations that the dimension converges to 3, with less assumptions. The situation in 4D 
may well be similar. 

It is interesting to note the possibility of measuring an effective average scalar curvature in 
a similar way. From eqn. ()2.11|) we see that, classically, the gradient of the spectral dimension 
function as a — > is related to this observable. We could define an effective scaling version of 
this observable, proportional to the gradient dD s (a)/da, and proceed in roughly the same way as 
in the AJL method of measuring the large scale dimension. For this to give a meaningful value, 
the range of values of a at which only the first two terms of the heat trace (|2.6p are contributing 
to the spectral dimension must be appreciably longer than the range over which quantum effects 
are important. From figure [2] we can see than our simulations do not allow this at present. 
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This would apply also to other ways of measuring curvture [251 [7] . This provides an interesting 
example of a situation in which concepts from general relativity and quantum field theory come 
together: although the average scalar curvature is defined as an integral of a function defined at 
each point, we know that this "maximally fine-grained" notion will no longer be the physically 
appropriate one in the quantum theory, and that we should only expect the concept of average 
scalar curvature to be useful in some coarse-grained sense, at some intermediate scales. In any 
case, the fit of the spectral dimension function to the scaled sphere shown in figure [2] gives good 
evidence for consistency with de Sitter beyond just the dimension or average scalar curvature. 

The two assumptions necessary to estimate the large scale dimension are not necessary to 
estimate the small scale dimension. The N = 200/c simulations are the largest available, and give 
a value consistent with 2 here, but inconsistent with 1.5. We will comment extensively on this in 
the following subsection. 

5.4 Discussion of the a — > limit 

It is interesting to compare the short-scale behaviour found here for the spectral dimension, to 
that derived from other approaches to quantum gravity. In most of the cases we are unable to 
give the analytic dependence on the diffusion time, and hence we can only discuss the qualitative 
behaviour, and speculate about the limiting value D s (a —> 0). While it would certainly be a bold 
extrapolation to claim that two approaches are describing the same physics just because they 
give the same number as a — > 0, it is intriguing to observe a certain universality of the results 
like it seems to happen in d = 4 [20]. In this respect it is useful to carry out this type of study 
in general dimensions, as for example d = 4 might be a special case in which several approaches 
agree on D s (a — > 0) despite having some fundamental difference. In this sense we think that 
a discrepancy in the results might have a stronger meaning, as a disagreement would probably 
signal some fundamental difference, and hence the dimension-dependence of the results might 
teach us something important. Our interest in three-dimensional CDT originates mainly from 
these concerns. 

We notice first of all that our result agrees with that found for d = 4 in [8], i.e. we find 
D s (a — > 0) ~ 2 both in three and four dimensions. We should have a word of caution here, in 
light of the recent observation [2] that, at least in the 4D case, the Planck length is of order 
the lattice spacing a at the currently possible parameter settings in Monte Carlo simulations, 
and remains fixed when the "bare Newton coupling" is fixed. This is something that we have not 
tested in our 3D simulations, but which could well hold true in this case. This observation implies 
that we are at the moment unable to probe physics well beyond the Planck scale in a reliable 
way (i.e. without discretization artifacts), hence the "<r —> 0" limit should be seen as a simple 
extrapolation of the results at a > L 2 p . In order to test whether such extrapolation is correct 
we would need to make the lattice length much smaller in comparison with the Planck scale, as 
described in [2] for 4D. This will not be attempted here, however. For the sake of discussion we 
will assume that we are seeing enough of the genuine UV behaviour at the values chosen for the 
parameters, and proceed to compare our CDT result to that coming from other approaches to 
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3D gravity. Further study will confirm whether this assumption is justified. 

In p3] the general a — > limit of the spectral dimension for asymptotically safe gravity in 
d dimensions was derived, resulting in D s {a — > 0) = d/2. Interestingly, for d = 2 this coincides 
with the AJL results, but for d = 3 it disagrees with the results found here. It should be stressed 
that in CDT, in both cases, the dimension is evaluated numerically, hence we can't be sure it 
will be exactly an integer. On the other hand the small errors give us confidence that a value of 
D s {a — > 0) = 3/2 in d = 3 can be excluded. 

In (n + l)-dimensional Hofava-Lifshitz gravity with characteristic exponent z one finds [19] 
that D s (a — > 0) = 1 + n/z. Hence D s (a — > 0) = 2 in the (2 + l)-dimensional case[f| with -2 = 2, 
which is the (2 + l)-dimensional analogue of the (3 + l)-dimensional case with z = 3 proposed in 
[26], n = z being the critical dimension of models characterized by z. This coincidence, together 
with the observation made in App. [Bj suggests that the link between CDT and Hofava-Lifshitz 
gravity deserves to be explored further. 

In the context of spin foam models of three-dimensional quantum gravity, the results of |18j 
give D s (a — > 0) = 2, but only after a transition at D s ~ 1.5 for small positive a. This behaviour 
could either be an artifact of the method used to determine the spectral dimension or something 
characteristic of spin foams. In any case, keeping in mind the existence of a minimal length in 
spin foam models, the limit a — > should probably be interpreted with some care. 

Finally, it is easy to adapt to d = 3 the calculation of [IS] for K-Minkowski and see that 
also in such case D s {a — > 0) = 2. As the dual of K-Poincare algebra, this specific type of non- 
commutative geometry has been suggested to play a role in quantum gravity [27]. We note that, 
as far as short scale spectral dimension is concerned, K-Minkowski agrees with CDT in d = 3 but 
disagrees in d = 4, a situation which seems to be opposite to that found about asymptotically 
safe gravity. 

In conclusion we see that studying the spectral properties at different spacetime dimensions 
can lead to some surprises and provide clues on potential relations between different approaches. 

6 Conclusions 

In this paper we have presented new results about the spectral dimension in CDT. On one hand 
we have applied to the three-dimensional case a short-scale analysis analogous to that of AJL, 
and on the other we have extended the analysis to larger scales. Both parts of this work represent 
an important extension to the understanding of the physics of CDT models. 

The short-scale results provide a hint and a first step towards an understanding of the in- 
dependence of the UV behaviour of CDT models. As explained in Sec. 15.41 more and more results 
are being produced about the spectral dimension in quantum gravity, and a comparison among 
them (as functions of the dimension d) can be very useful. 

The large-scale results represent striking evidence for classical behaviour at large scales, which 
goes beyond the time-dependence of the spatial volume [6], as the spectrum of the Laplacian 

5 Note that n = 2 here corresponds to what we have called d = 3 in the rest of the paper. 
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contains far more detailed information than that. The least that can be said is that evidence 
presented here indicates that in the N — > oo limit (at fixed kq) the spectral dimension function 
tends to that of a stretched sphere. By the freedom we have to rescale the "cosmological" proper 
time, we conclude that results are consistent with the production of a de Sitter ground state in the 
model. We note again that this spectral dimension function contains a great deal of geometrical 
information, which can completely characterise some simple geometries. Ideally, we would like 
to say that, once the CDTs are appropriately scaled, the spectral dimension function tends to 
that of a sphere with standard metric. This would provide evidence for a necessary and sufficient 
condition for the ground state of three-dimensional CDTs to be exactly a sphere in the classical 
limit. We can establish this using the 1-to-l relation between the spectral dimension function 
and spectrum of the Laplacian, and the well-known theorems [TU] mentioned in section [2] (at least 
up to some remaining technical concerns involving smoothness assumptions). The comparison to 
the stretched sphere prevents such a clean statement, but still gives strong evidence in favour the 
same conclusion. Further, in the real world there are quantum corrections at short scales that 
we do not wish to remove by such a limiting procedure. It is tempting to conjecture that, if the 
spectral dimension function of a geometry approximates to that of a sphere at large scales, the 
the geometry is "approximately isometric to the sphere at large scales" in some appropriate sense; 
this would again allow a cleaner statement. But even without such an ambitious conjecture, the 
match of the spectral dimension functions exhibited in figure [2] presents powerful new evidence 
for the approximation of 3D CDT simulation results to de Sitter space at large scales. 

Extensions of the work could involve more extensive comparisons with the spectral dimension 
functions arising from other approaches to quantum gravity. This kind of comparison could 
suggest ways to derive an estimation of the Planck length (or some proportional scale) from the 
small scale form of the spectral dimension function for CDTs, for example by relating the Planck 
length to the location of the maximum of D s (a). This would be useful if and when studies 
are undertaken to investigate the sub-Planckian regime of the models. Such studies would be 
necessary in order to assess the existence of a non-trivial continuum limit, i.e. one with local 
degrees of freedom and a finite Planck length. In addition, it would be interesting to see if such 
limit could be taken in a way compatible with the rescaling needed to have an isotropic universe: 
this would shed light on the connections between CDTs and Hof ava-Lifshitz gravity. 

In view of the successful application of our method to 3D CDT, the same techniques can now be 
applied to the more physically interesting 4D CDT models. There are no new conceptual problems 
to overcome to do so; the extension merely involves an application of the same procedures to the 
previously studied 4D simulations. This will provide interesting new ways to explore the large 
scale geometry, as well as to investigate small-scale behaviour as we push simulations further into 
the sub-Planckian regime. 
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Appendices 

A Stretched sphere 

A.l Definition and geometric properties 

We define a three-dimensional stretched sphere Sf by uniformly stretching the proper distance 
in the longitudinal direction. Its metric is 

ds 2 s s = r 2 (s 2 dip 2 + sin 2 ip (d9 2 + sin 2 9 d<p 2 )) , (A.l) 

where ip,9 G [0,7r] and cp G [0, 2tt], and s is the deformation parameter (for s = 1 we have the 
standard metric on a sphere of radius r). It is easy to see that a manifold with such metric has 
conical singularities at the poles; for example, we can approximate s'mip ~ ip near the North pole, 
and by making the substitution {p = rstp, V 9 = (also restricting to the subspace 9 = ir/2 for 

simplicity) , we find the metric of a cone 



ds 2 = dp 2 + p 2 d<p 2 , (A.2) 



where ip 6 [0, 2ir/s]. 

The scalar Ricci curvature is 



6sinV + 2(s 2 -l) 

K = 22 -2, ' \ A - 6 ) 

r z s A sin ip 



and is divergent as expected in ip = 0, ir (incidentally its integral is finite, but that is not true for 
higher order invariants). On a manifold with conical singularities no closed analytical expression 
is known for the coefficients of the heat kernel expansion, but we will be able to compute directly 
the spectrum for this particular case. 
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An embedding in 



can be obtained with the coordinate parameterization 
x o = r y s 2 — cos 2 i/j' dip' , 



x\ = r cos (p sin sin ip , 
X2 = r sin cfi sin 9 sin ip , 
X3 = r cos # sin ^ . 

A section of such embedding is shown in FigJBJ 



(A.4) 




Figure 6: An embedding in R 3 of the 8 = ir/2 section of the three-dimensional stretched sphere, with 
r = 1.2 and s = 1.96. 



A. 2 The Laplacian and its spectrum 

The scalar Laplacian associated to the metric (jA.ip is 

A S'if = 2 2 1 2 > ( sin2 V> <W) + 2 ■ 2~i ■ fl d e (sin6> <%/) + . * . 2 
a r z s z sm ^ rsin ipsmU r z sm ipsm V 



1 , 1 ( A - 5 > 

= 2 2 -2 7 ^(sm 2 ^^/)+ . 2 - Affl/, 

r z s z sm ^ r z sin ?/> 

where A S 2 is the Laplacian on a unit 2-sphere. In order to look for its spectrum we start with a 
separation of variables, i.e. we make the ansatz 

f(ip,e,<j>) = P (ip)Yr(e,cf>), (A.6) 

where Yj m {9, <p) are spherical harmonics (I = 0, 1, 2... and m = —I, —I + 1, ...I), so that 



A s ,Yr(e, </>) = (a 2 + c ^d e + -J^if^, 0) = -«(« + i)*r(*, 0) • (A.7) 

sm f sin w 
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As a consequence, the eigenvalue equation reduces to an ordinary differential equation for p(ijj): 

p{i>) = -E p{i>) . (A.8) 



1 



1 2 2cosV„ 1(1 + 1) 

9 U tt> • 9 • / U ll> ■ 1 , 



(A.10) 



The generic solution for this equation is found to be of the form 

VSirr^ 

where Pjr and are the associated Legendre functions of the first and second kind, with A = 
2Vi+g^r 2 _E zl i an( ^ ^ _ ly/J _|_ 4s 2 ^/ _|_ ]^ anc j coefficients ^4 and i? to be fixed by boundary 

conditions. Imposing regularity at ip = we find that 

A 7T COS fill 

B 2 sin fiir ' 

while B can be fixed by normalization. 

Finally, imposing regularity also at ip = tt, the eigenvalues E are restricted to 

2j(j + 1 - 2/ + 2/x) + 2Z(Z - 1 + s 2 (/ + l) - 2 /i) + 2/i - 1 

% = ^2-2 > (A" 11 ) 

with j > I a non-negative integer. The degeneracy of the eigenvalues is that inherited from the 
spectrum over S 2 , i.e. 

Dji = 2l + 1. (A. 12) 

For s = 1 the eigenvalues become independent of I and the spectrum of the round 3-sphere is 
recovered: 

f:r i(i+2> - .A-, 

Once we have the full spectrum the heat kernel trace can be directly evaluated by the formula 

^w = ^EE^'. ( A - 14 ) 

s s j=o 1=0 

with the volume given by Vg3 = 2ir 2 r s s. 



B Solution to (2+l)-dimensional z = 2 Horava-Lifshitz gravity 

A recent proposal by Petr Hofava |26j of a power-counting renormalizable theory of gravity in 
3+1 dimensions has attracted lot of attention by the physics community. A possible link between 
this model and CDT was pointed out in |19j . and we wish to explore it further here in light of 
our results. In general, in n + 1 dimensions, the idea is to construct an anisotropic theory of 
gravity, characterized by a critical exponent z = n, such that power-counting renormalizability is 
guaranteed without introducing higher-derivatives in time (which would typically spoil unitarity). 
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The case of our interest is hence a z = 2 theory in 2+1 dimensions, which has already been 
considered in [28]. 

We will consider here a version without detailed-balance condition and in Euclidean signature, 
for which the action is 

S = j dtd 2 xN^g j-^ XK 2 - KijK ij - 2A + - 7^ j , (B.l) 

where g is the determinant of the spatial metrical, ( 2 lR its Ricci scalar, N the lapse function, 
Kij the extrinsic curvature of the leaves of the foliation, and K its trace. The coupling k 2 is 
proportional to Newton's constant, and A is the cosmological constant, while A and 7 characterize 
the deviations from full diffeomorphism invariance (1 = 1 and 7 = corresponding to general 
relativity in 2+1 dimensions). In the following, since we are interested in the infrared limit of 
the theory, we will restrict to 7 = 0, hence leaving only to A / 1 the characterization of the 
anisotropy. 

The general equations of motion can be easily found and are similar to those in [29], but we 
will not report them here. We are interested in spatially spherically-symmetric solutions, and 
therefore we make the ansatz 

ds 2 = N 2 {il>) dV 2 + fW(d9 2 + sin 2 9 d(f> 2 ) . (B.2) 



As a consequence the equations of motion reduce to 



2X ~ 1 ft) 2 + 2A--=0, (B.3) 



2iV 2 \fj f 
2A-1 (I ff'V N' f f" 



, , , + — — -— -AAT = 0. (B.4) 
2N \2\f J N f f ) 1 > 

It is easily verified that a solution is given by 

/0/0 = ]-sin 2 ^, (B.5) 

N 2 W = ^-±. (B.6) 

This corresponds exactly to the metric (|A. 1|) for r 2 = 1/A and s 2 = 2A — 1. 

We hence see that the large scale spectrum found in our CDT simulations in general matches 
the above solution to Hofava-Lifshitz gravity. As explained in the main text we are at this stage 
free to translate our results to the continuum language in such a way that the round sphere is 
matched, but we should keep in mind the general result when taking into account the scaling 
of kq in future studies: for example, when taking the continuum limit it might turn out to be 
impossible to keep the Planck length finite and at the same time have a round sphere at large 
scales. In the renormalization group language this would mean that the coupling A does not flow 
to A = 1 in the infra-red. We hope to come back to this issue in the near future. 



^Despite being in Euclidean signature we keep using a spacetime terminology, for lack of a better one. 



23 



References 



S. W. Hawking, The path integral approach to quantum gravity. In 'General Relativity' ed. 
S.W. Hawking, W. Israel, Cambridge Univ. Pr. 

[2] J. Ambjorn, J. Jurkiewicz and R. Loll, Quantum gravity as sum over spacetimes, 
larXiv: 0906. 39471 

[3] J. Ambj0rn, J. Jurkiewicz and R. Loll, Reconstructing the universe, Phys. Rev. D72 (2005) 
064014 |arXiv : hep-th/0505 154] . 

[4] J. Ambjorn, J. Jurkiewicz and R. Loll, Emergence of a J^D world from causal quantum 
gravity, Phys. Rev. Lett. 93 (2004) 131301 [ arXiv : hep-th/0404156) . 

[5] J. Ambj0rn, B. Durhuus and T. Jonsson, Quantum geometry. A statistical field theory 
approach. Cambridge Univ. Pr., 1997. 

[6] J. Ambjorn, A. Gorlich, J. Jurkiewicz and R. Loll, The Nonperturbative Quantum de Sitter 
Universe, Phys. Rev. D78 (2008) 063544 [arXiv: 0807. 448 lj . 

[7] J. Henson, Coarse graining dynamical triangulations: a new scheme, Class. Quant. Grav. 
26 (2009) 175019 [arXiv: 0907 .5602] . 

[8] J. Ambj0rn, J. Jurkiewicz and R. Loll, Spectral dimension of the universe, Phys. Rev. Lett. 



95 (2005) 171301 |arXiv: hep-th/0505 113 



[9] D. V. Vassilevich, Heat kernel expansion: User's manual, Phys. Rept. 388 (2003) 279-360 
|arXiv : hep-th/0306 138| . 

[10] S. Tanno, Eigenvalues of the Laplacian of riemannian manifolds, Tohoku Math. Journ. 25 
(1973) 391-403. 

[11] S. Weinberg, Ultraviolet Divergences in Quantum Theories of Gravitation. In 'General 
Relativity' ed. S.W. Hawking, W. Israel, Cambridge Univ. Pr. 

[12] M. Niedermaier and M. Reuter, The Asymptotic Safety Scenario in Quantum Gravity, 
Living Rev. Rel. 9 (2006) 5. 

[13] R. Percacci, Asymptotic Safety. In 'Approaches to Quantum Gravity: Towards a New 
Understanding of Space, Time and Matter' ed. D. Oriti, Cambridge Univ. Pr. 
larXiv: 0709 .38511 

[14] O. Lauscher and M. Reuter, Fractal spacetime structure in asymptotically safe gravity, 
JHEP 10 (2005) 050 |arXiv : hep-th/0508202] . 

[15] D. Benedetti, Fractal properties of quantum spacetime, Phys. Rev. Lett. 102 (2009) 111303 
[arXiv:0811. 



24 



[16] L. Modesto, Fractal Structure of Loop Quantum Gravity, arXiv : 0812 . 2214, 
[17] L. Modesto, Fractal Quantum Space-Time, arXiv : 0905 . 1665, 

[18] F. Caravelli and L. Modesto, Fractal Dimension in 3d Spin-Foams, arXiv : 0905 . 2170, 

[19] P. Horava, Spectral Dimension of the Universe in Quantum Gravity at a Lifshitz Point, 
Phys. Rev. Lett. 102 (2009) 161301 [arXiv: 0902. 3657] . 

[20] S. Carlip, Spontaneous Dimensional Reduction in Short- Distance Quantum Gravity?, 
larXiv : 0909 . 33291 

[21] J. Ambj0rn, J. Jurkiewicz and R. Loll, Dynamically triangulating Lorentzian quantum 
gravity, Nucl. Phys. B610 (2001) 347-382 |arXiv : hep-th/0 105267] . 

[22] J. Ambjorn, J. Jurkiewicz, R. Loll and G. Vernizzi, Lorentzian 3d gravity with wormholes 
via matrix models, JHEP 09 (2001) 022 | arXiv : hep-th/0106082 . 

[23] D. Benedetti, R. Loll and F. Zamponi, (2+1) -Dimensional Quantum Gravity as the 
Continuum Limit of Causal Dynamical Triangulations, Phys. Rev. D76 (2007) 104022 
|arXiv: 0704 .3214] . 

[24] J. Ambj0rn, J. Jurkiewicz and R. Loll, Non-perturbative 3d Lorentzian quantum gravity, 
Phys. Rev. D64 (2001) 044011 |arXiv : hep-th/001 1276] . 

[25] B. V. de Bakker and J. Smit, Curvature and scaling in 4-d dynamical triangulation, Nucl. 
Phys. B 439 (1995) 239 |arXiv : hep-lat /94070I4] . 

[26] P. Horava, Quantum Gravity at a Lifshitz Point, Phys. Rev. D79 (2009) 084008 
[arXiv: 090 1.3775] . 

[27] G. Amelino-Camelia, L. Smolin and A. Starodubtsev, Quantum symmetry, the cosmological 
constant and Planck scale phenomenology, Class. Quant. Grav. 21 (2004) 3095-3110 
|arXiv : hep-th/0306134] . 

[28] P. Horava, Membranes at Quantum Criticality, JHEP 03 (2009) 020 |arXiv: 0812. 4287] . 

[29] H. Lu, J. Mei and C. N. Pope, Solutions to Horava Gravity, Phys. Rev. Lett. 103 (2009) 
091301 [arXiv: 0904^1595] . 



25 



